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INTRODUCTION 


Military  requirements  for  more  efficient  smoke  munitions  have  led  to  the 
development  of  the  XM-825  felt  wedge  white  phosphorus  (FWWP)  155-mm  projec¬ 
tile.  The  basic  goals  of  the  XM-825  were  to  provide  a  smoke  munition  that 
nearly  eliminated  the  exothermic  and  the  quasi -instantaneous  source  character¬ 
istics  associated  with  bulk-filled  white  phosphorus  (WP).  The  XM-825  projec¬ 
tile  accomplishes  this  by  using  92  phosphorous-impregnated  felt  wedges*  that 
provide  an  area  smoke  source  of  a  continuous  nature  with  little  or  no  exother¬ 
mic  effect  when  dispersed. 

The  continuous  area  source  has  introduced  some  complexities  to  the  diffusion 
modeling  associated  with  this  candidate  munition.  However,  several  area 
source  approaches  to  the  downwind  dispersion  of  windborne  material  are  avail¬ 
able.  Notable  among  these  hypotheses  are  a  pseudopoint  model  advocated  by 
Turner,1  an  infinite  line  solution  also  suggested  by  Turner,  and  the  algebraic 
integration  method  suggested  by  Pasquill.2 

The  three  approaches  to  area  source  diffusion  were  investigated  within  the 
framework  of  the  KWIK  smoke  algorithm.3  Results  obtained  during  evaluation  of 
the  three  methods  using  data  from  Dugway  Proving  Ground4  were  encouraging. 


THEORETICAL  CONSIDERATIONS 


The  XM-825  Candidate  Smoke  Munition 

The  XM-825  is  a  155-mm  projectile  designed  to  disperse  92  phosphorous- 
saturated  felt  wedges*  which  will  produce  a  screening  smoke  for  5  to  10  min. 
The  total  warhead  payload  is  7.48  kg,  consisting  of  1.13  kg  of  felt  wedges  and 
6.35  kg  (+0.10  kg)  of  white  phosphorus.  Preliminary  estimates  indicate  that 
this  configuration  is  82  percent  efficient  by  weight  and  produces  effective 
obscuration  in  the  visible  portion  of  the  electromagnetic  spectrum  for  up  to 
450  s. 


♦Configuration  of  XM-825  as  of  March  1980. 

burner,  D.  B.,  Workbook  of  Atmospheric  Dispersion  Estimates,  Public  Health 
Service  Publication  No.  999-AP-2 6,  Cincinnati,  Ohio,  Rev.  1969. 

2Pasquill,  F.,  A tmospheri c  D if fusion ,  second  edition,  Halsted  Press  Div,  John 
Wiley  and  Sons,  Inc.,  New  Vork,  429  pp,  1974. 

3Umstead,  R.  K.,  R.  Pena,  and  F.  V.  Hansen,  KWIK:  An  Algorithm  for 

Calculating  Munition  Expenditures  for  Smoke  Screening/Observation  in  Tactical 
Situations,  ASL-TR-O03O,  19)9. 

40ugway  Proving  Ground,  Development  Test  1  of  155  irrn  Smoke  Projectiles  (XM-803 
and  XM-825),  Defense  Technical  Information  Center,  Cameron  Station, 
Alexandria,  Virginia,  1979. 
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The  Pseudopoint  Solution 


One  approach  to  solving  area  source  diffusion  problems  is  to  assume  that  the 
area  can  be  represented  by  a  pseudopoint  located  at  some  distance  xQ  upwind  of 
the  center  of  the  area.  Thus,  the  origin  of  the  vertical  dispersion  parameter 
o2  is  at  the  pseudopoint--in  this  case,  the  upwind  edge  of  the  XM-825  burst 
pattern.  The  vertical  dispersion  parameter  (standard  deviation)  can  then  be 
written  as: 


o 


z 


(1) 


where  xQ  is  the  burst  radius  and  x  the  downwind  distance.  The  coefficients 
and  indices  c  and  d  for  equation  (1)  are  dependent  upon  both  surface  roughness 
and  the  stability  of  the  surface  boundaty  layer  according  to  Pasquill.*  The 
obscuration  power  of  a  screening  smoke  is  dependent  upon  the  crosswind  inte¬ 
grated  concentration  (CVI1C)  of  the  aerosol  in  an  optical  path  traversing  the 
plume.  The  integrated  concentration  through  a  nonbuoyant  continuous  plume  can 
be  shown  to  be  given  by: 


where  x  is  a  munition  efficiency,  ft  a  relative-humidity-dependent  yield 
factor,  Q  the  source  strength  in  grams  per  second,  Y  the  mean  windspeed 
in  meters  per  second,  and  z  the  height  above  the  surface  of  the  optical 
path.  Equation  (2)  is  an  empirically  modified  version  of  a  Gaussian  plume  in 
that  a  power  of  3/2  is  used  in  place  of  the  more  conventional  second  power. 
This  is  done  to  compensate  for  nori-Gaussian  effects  upon  vertical  dispersion, 
based  upon  experimental  results  obtained  by  Pasquill,2  Huang,5  and  Horst.6 


2Pasquill,  F.,  Atmospheric  Diffusion,  second  edition,  Halsted  Press  Div,  John 
Wiley  and  Sons,  Inc.,  New  York,  429  pp,  1974. 

sHuang,  C.  H.,  "A  Theory  of  Dispersion  in  Turbulent  Shear  Flow,"  Atmospheric 
Environ,  13:453,  1979. 

6Horst,  T.  W.,  "Lagrangian  Similarity  Modeling  of  Vertical  Diffusion  From  a 
Ground-Level  Source,"  J  Appl  Meteorol,  18:733,  1979. 
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The  Infinite  Line  Source 


Turner1  suggests  that  area  sources  may  be  treated  by  considering  the  area  as  a 
line  source  by  compressing  all  of  the  sources  into  an  infinite  line  running 
through  the  axis  of  the  burst  pattern.  Concentrations  at  a  given  distance 
downwind  may  then  be  calculated  from: 


x  = 


sin 


2qA  a 


<J>  /Z7  o^V 


—  exp 


(3) 


where  <j>  is  the  angle  between  the  wind  direction  and  the  line  source,  and  q  = 
Q/L  with  L  being  the  diameter  of  the  burst.  Equation  (3)  is  not  valid  when  *> 
is  less  than  45°.  The  vertical  standard  deviation  az  is  calculated  in  the 
same  fashion  as  in  the  pseudopoint  model.  The  horizontal  standard  deviation 
Oy  is  required  to  calculate  CVIIC  and  is  determined  from: 


ax 


(4) 


where  the  coefficient  a  and  index  b  are  functions  of  the  Pasquill  Stability 
Category.  CWIC  can  then  be  calculated  from  the  following: 


xcwic  =  X^L  +  4'3  °y)  (5) 

The  factor  4.3  is  included  to  delineate  the  visible  edge  of  the  plume,  where 
the  concentration  decreases  to  10  percent  of  the  centerline  value. 

A  comparison  was  made  between  the  infinite  line  source  and  a  finite  line 
source,  which  is  basically  the  same  equation  with  an  edge  effect  error  func¬ 
tion.  The  difference  in  downwind  concentrations  at  screening  distances  were 
found  to  be  insignificant. 


Algebraic  Integration  Method 

The  pseudopoint  and  infinite  line  source  models  are  very  simplistic  models 
based  on  unrealistic  geometries  of  the  burst  pattern.  They  are  only  capable 
of  predicting  an  average  value  for  CWIC.  A  more  sophisticated  approach  would 
be  to  model  the  geometry  of  the  burst  to  give  a  more  realistic  answer.  The 
algebraic  integration  approach  of  Pasquill2  attempts  to  model  the  geometry 
with  a  circular  burst  pattern,  as  shown  in  figure  1.  The  circular  pattern  is 


burner,  D.  B.,  Workbook  of  Atmospheric  Dispersion  Estimates,  Public  Health 
Service  Publication  No.  999-AP-26,  Cincinnati,  Ohio,  Rev.  i9b9. 

2Pasquill,  F.,  Atmospheric  Diffusion,  second  edition,  Halsted  Press  Div,  John 
Wiley  and  Sons,  Inc.,  Mew  York,  429  pp,  1974. 
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divided  into  crosswind  strips  of  equal  downwind  lengths  ax.  The  strips  are  of 
varying  widths,  with  the  submunitions  of  the  XM-825  spread  uniformly  through¬ 
out  the  pattern.  A  weighting  factor  W  is  used  to  represent  the  relative 
importance  of  the  separate  contributions  from  each  strip  and  is  given  by: 

W  .  nd-<0  .  („  .  !)(!-«)  ,6, 


where  n  is  the  number  of  the  strip  (the  first  strip  being  farthest  down¬ 
wind).  The  index  d  is  identical  with  that  of  equation  (1).  A  sample  of  the 
weights  for  six  strips  for  a  surface  roughness  length  of  10  cm  and  a  Pasquill 
Stability  Categoy  of  "B"  is  as  follows: 


n  =  1  2  3  4  5  6 

weighting  factor  =  1.0  0.11  0.07  0.05  0.04  0.04 


Eighty-five  percent  of  the  source  is  contained  in  the  first  two  strips  due  to 
the  rapid  decline  in  the  contribution  of  the  more  distant  strips. 

The  contribution  to  the  downwind  concentration  of  a  given  strip  n  is: 


X  ■  w  A  qn{l-a) 


-1 


inAx 


(n-l)Ax 


(7) 


where  A  =  -  ,  q  is  the  source  strength  of  strip  n  in  grams  per 

sin  $  it 

meter  per  second  and  j  is  a  dummy  variable.  Calculating  the  total  CWIC  at 
some  downwind  sampling  distance  requires  that  concentrations  first  be  calcu¬ 
lated  for  each  separate  strip  and  then  be  summed  over  the  six  strips,  using 
the  appropriate  weighting  factors.  In  addition,  oz  must  be  calculated  sepa¬ 
rately  for  each  strip,  since  the  successive  strips  are  farther  upwind  by  a 
distance  of  ax.  CWIC  is  then  evaluated  from: 


xcwic  =  x^Ls  +  4*3  °y)  (8) 


assuming  that  Ls  is  the  width  of  the  source  strips. 
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MODEL  RESULTS  AND  EVALUATIONS 

The  three  area  sources  postulates  were  used  to  evaluate  seven  sets  of  experi¬ 
mental  data  extracted  from  data  reports  published  by  Dugway  Proving  Ground.4 
Of  the  seven  examples,  three  were  chosen  for  inclusion  in  this  report.  The 
remainder  were  eliminated  because  of  apparent  improper  burning  characteris¬ 
tics,  unfavorable  wind  direction,  or  a  combination  of  both.  A  summary  of  test 
conditions  is  listed  in  table  1  and  the  test  grid  is  illustrated  in  figure  2. 

The  integrated  concentration  values  of  the  pseudopoint,  infinite  line,  and 
algebraic  integration  models  are  as  follows: 


Algebraic  Integration 


Test 

Pseudopoint 

Infinite  Line 

Peak 

Mean 

G1-A2 

0.66  g  m"  2 

1.62  g  m~2 

2.19  g  m'2 

1.40  g  m“2 

G2-A1 

0.50  g  m“2 

1.24  g  m-2 

1.99  g  in-2 

1.34  g  m-2 

G4-A1 

0.49  g  m~2 

1.20  g  m'2 

2.11  g  m-2 

1.38  g  m*2 

The  pseudopoint  model  generates  average  integrated  concentrations  which  are 
consistently  the  lowest  of  the  three  models.  It  underpredicts  the  measured 
values,  as  shown  in  figures  3,  4,  and  5.  This  is  not  surprising,  due  to  the 
simplifying  assumptions  in  the  model. 

The  infinite  line  model  does  better  than  the  pseudopoint  model  in  giving 
average  integrated  concentrations  that  are  closer  to  the  measured  values,  as 
shown  in  figures  3,  4,  and  5.  However,  during  the  final  150  s  of  burn  time 
the  model  overpredicts.  For  test  G4-A1,  the  results  are  slightly  high  for  the 
entire  burn. 

The  algebraic  integration  model  (which  will  be  called  the  circle  model  due  to 
the  shape  of  the  modeled  burst)  gives  the  best  results  of  the  three.  This 
model  generates  both  a  peak  and  a  geometric  mean  value  for  CWIC.  The  peak 
value  is  the  highest  CWIC  that  is  expected  for  the  given  meteorological 
conditions  and  the  distance  from  the  optical  line  of  sight.  The  threshold 
value  of  smoke  which  will  obscure  the  visible  wavelengths  is  calculated  from: 


CL 


fcn  T 


smk 


-a 


(9) 


4Dugway  Proving  Ground,  Development  Test  1  of  155  mm  Smoke  Projectiles  (XM-803 
and  XM-825),  Defense  Technical  Information  Center,  Cameron  Station, 
Alexandria,  Virginia,  1979. 
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where  CL  is  the  concentration,  Tsmk  the  transmittance,  and  a  the  extinction 
coefficient.  The  threshold  concentration  is  calculated  to  be  0.7  g  m'2.  This 
value  occurs  at  an  average  burn  time  of  450  $.  By  assuming  the  peak  value 
occurs  at  the  beginning  of  the  burn,  an  estimate  of  the  emission  rate  curve  of 
the  XM-825  can  be  made.  The  results  are  labeled  "circle"  in  figures  3,  4,  and 
5.  The  actual  emission  characteristics  curve,  as  reported  by  Carter/  indi¬ 
cates  that  the  emission  rate  of  the  XM-825  decreases  with  time  after  its 
initial  buildup  time.  This  is  due  to  the  peculiarities  of  the  bum  charac¬ 
teristics  of  the  phosphorous-impregnated  felt  submunitions.  (The  geometric 
mean  value  is  listed  with  the  model  results  shown  previously.) 

The  overestimation  of  CWIC  by  this  model  for  the  first  75  s  of  the  burn--which 
is  evident  in  figures  3,  4,  and  5 — is  a  function  of  the  buildup  time,  the  time 
required  for  the  smoke  plume  to  reach  the  measuring  line,  and  the  percentage 
of  the  cloud  that  actually  drifts  through  the  optical  line  of  sight.  As  with 
the  infinite  line  model  there  is  a  slight  overestimation  of  CWIC  in  test 
G4-A1.  Both  of  these  overestimates  could  be  caused  by  the  failure  of  some  of 
the  submunitions  to  burn  properly,  thus  giving  lower  concentrations  than  would 
normally  be  expected. 

The  preliminary  center-to-center  impact  separations  can  be  computed  based  on 
the  weighted  circular  model.  The  equation  to  compute  the  separation  is  a 
variation  of  equation  (7),  which  is  as  follows: 


x 


n 

zfw  q 

L  n  " 

,v  c  xcwic 


(10) 


For  stability  category  C  the  impact  separation  is  132  m,  and  for  category  B 
the  impact  separation  is  121  m. 


CONCLUSIONS 

The  pseudopoint  model  underpredicts  the  average  integrated  concentrations  in 
all  cases  investigated.  Therefore,  this  model  would  overpredict  the  number  of 
munitions  necessary  to  screen  and  would  not  be  ah  acceptable  model.  The 
infinite  line  model,  although  better  than  the  pseudopoint  model  still  over¬ 
predicts  the  final  150  s  of  the  burn.  This  would  cause  too  few  munitions  to 
be  expended,  risking  the  deterioration  of  the  screen.  The  circle  model 
appears  to  be  the  best  fit  to  the  experimental  data.  It  gives  a  fairly  good 


7Dugway  Proving  Ground,  Basic  Smoke  Characterization  Test,  TEC0M-DPG-TP-77- 
311,  Dugway  Proving  Ground,  Utah,  1975. 
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estimate  of  the  emission  curve,  as  shown  in  figures  3,  4,  and  5.  The  over¬ 
estimation  of  the  concentrations  in  the  model  only  occurs  during  the  first 
75  s  of  the  burn,  which  is  the  buildup  time  of  the  smoke  generated  by  this 
munition.  This  model  has  the  promise  of  giving  good  expenditure  estimates, 
and  it  deserves  further  investigation. 


■ 


TABLE  1.  SUMMARY  OF  TEST  CONDITIONS 


Test 

G1-A2 

G2-A1 

Date  -  Time 

8  Aug  78  -  1020 

8  Aug  78  -  1 

Wind  direction  (8  m) 

339° 

338° 

Windspeed  (8  m) 

3.8  m/s 

3.8  m/s 

Relative  humidity  (4  m) 

16% 

14% 

Temperature  (8  m) 

29.3°C 

30. 8° C 

Visibility  (km) 

48 

48 

Cloud  cover 

0 % 

0% 

Stability  category 

C 

B 

Entry  angle 

15.9° 

15.9° 

Burst  altitude  (m) 

102 

103 

Time  to  obscuration  (s) 

9.5 

11 

12 


A  MID  1R  TRANSMITTERS 
•  MID  IR  RECEIVERS 
O  GRID  CENTER 

Figure  2.  Test  grid  for  XM-825  FWWP  snoke  mini tl on 


CLQf  m  2  ^ 


Figure  3.  CHIC  for  XM-825  FHHP  smoke  munition.  Test  G1-A2. 
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